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ABSTEACT 


This  report  proposes  a  modification  to  the  Wheeler  fatigue  crack 
retardation  model.  The  modified  model  is  derived  and  then  used  to  predict 
existing  data  for  the  number  of  delay  cycles  following  a  single  overload 
in  2024— T3  Aluminum,  4340  Steel,  and  Tl— 6A1-4V  Titanium,  alloy.  All 
predictions  were  within  essentially  a  factor  of  two  of  the  experimental 
data. 
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SYMBOLS 


Aa  -  current  crack  growth  increment  since  overload 

C,  n  -  crack  growth  ,  rate  constants 

da/dN  -  crack  growth  rate 

^max  ~  maximum  stress  intensity 

K^in  ”  minimum  stress  intensity 

Kql  ~  overload  stress  intensity 

AK  -  stress  intensity  range  (K  -  K  .  ) 

max  min 

AK^££  -  effective  stress  intensity  range  . 

AK^I^  -  threshold  stress  intensity  range 

m  -  Wheeler  shaping  exponent 

-  observed  delay  cycles 

Np  r-  predicted  delay  cycles 

R  -  K  .  /K 

min  max 

S  -  crack  arrest  factor 

-  overload  interaction  zone 

^ys  -  yield  stress 
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I  INTRODUCTION 

The  qualitative  effect  of  an  overload  on  a  fatigue  crack  propagating 
at  some  lower,  cyclic  load  is  well  known.  A  tensile  overload  will  delay 
or  retard  subsequent  fatigue  crack  growth  below  that  expected  for  the  case 
of  no  overload.  Any  crack  growth  analysis  that  does  not  accomt  for  this 
load  interaction  will  predict  an  overly  conservative  crack  growth  life. 

The  Wheeler  model  was  one  of  the.  first  attempts  to  quantify  fatigue 
crack  retardation  (Ref  1).  The  biggest  drawback  to  this  model  is  the 
empirical  ’’shaping  exponent”,  m.  Forcing  the  exponent  m  to  be  a  constant 
may  provide  a  good  correlation  between  the  model  and  a  particular  set  of 
experimental  data,  but  an  important  loss  of  generality  is  incurred.  This 
report  proposes  a  modification  to  the  Wheeler  model  which  allows  the  model 
to  be  used  without  reliance  on  data  fitting  and  without  the  subsequent 
limiting  to  a  specific  set  of  loading  parameters. 


II  MODEL  DEVELOPMENT 

Gallagher  (Ref  2)  presented  the  Wheeler  model  in  the  following  stress 
intensity  format : 
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Experiments  involving  single  overloads  in  2024-T3  Aliminum  (Ref  3) 
and  Ti-6A1-4V  Titanium  alloy  (Ref  4)  indicate  that  there  is  a  particular 


value,  S,  of  K„, /K  such  that  when  K-., /K  -  S,  crack  arrest  occurs. 
’  *  OL  max  OL  max 
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For  2024-T3  (and  possibly  other  altaminum  alloys),  S  is  about  2.3.  For 
Ti-6A1-4V,  S  is  about  2.8. 

The  limiting  condition  for  crack  arrest  is  /K  =  S.  Immediately 

OL  max 

after  the  overload,  Aa  =  0.  Substituting  these  relations  into  Eq.  3  yields 

2m 

AK^^^  (at  arrest)  =^l|  (4) 

Gallagher  (Ref  2)  stated  that,  for  the  limiting  case  of  crack  arrest, 
the  effective  stress  intensity  range  equals  the  threshold  stress  intensity 
range.  Using  this  concept  and  the  modified  Wheeler  formulation  for 
(at  arrest)  [Eq  (4) ]  leads  to 

2m 

hi  " 

-[sj  [“] 

Solving  Eq  (5)  for  the  exponent  m. 


Thus,  it  is  evident  that  the  Wheeler  exponent  m  is  not  a  constant  but 
depends  on  the  material  and  the  loading  subsequent  to  the  overload. 

Ill  DATA  CORRELATION 

The  number  of  delay  cycles  (N,)  following  a  single  overload  has  been 

a 

observed  for  2024-T3  (Ref  3,  5,  6),  4340  Steel  (Ref  2),  and  Ti-6A1-4V 
(Ref  4,  7).  A  delay  cycle  is  defined  as  any  post-overload  cycle  during 
which  the  growth  rate  is  less  than  that  expected  for  the  case  of  no  overload. 
See  Fig.  1.  In  the  2024-T3  and  4340  experiments,  load  shedding  techniques 
were  employed  to  obtain  constant  AK  loading  before  and  after  an  overload. 

See  Fig.  2.  Because  of  slightly  different  experimental  techniques,  the 
observed  delay  cycles  in  the  Ti-6A1-4V  experiments  correspond  only  approxi¬ 
mately  to  Nj  as  defined  in  Fig.  1. 
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Loading  parameters  and  observed  delay  cycles  from  Eeferences  2 
through  7  are  summarized  in  Table  I.  Eq  (1)  in  conjunction  with  Eq  (6) 
has  been  \ised  in  a  numerical  integration  routine  to  predict  the  number 
of  observed  delay  cycles  for  the  various  sets  of  loading  parameters  in 
Table  I.  Fig.  3  shows  the  correlation  between  predictions  (N  )  and  actual 

p 

data  (N^)« 

The  arrest  factors  S  for  2024-T3,  4340,  and  Ti-6A1-4V  were  assumed 
to  be  2.3,  2.3,  and  2.8,  respectively.  Table  IT.  summarizes  the  growth 
rate  constants,  C  and  n,  that  were  used  for  the  various  materials  and  R 
ratios.  These  constants  were  determined  by  a. least  squares  fit  to  da/dN 
versus  AK  data  presented  in  References  2  through  7,  or,  where  applicable, 
the  constants  were  taken  directly  from  these  references.  Threshold  stress 
intensity  ranges  for  2024-T3,  4340,  and  Ti-6A1-4V  were  assumed  to  be 
2,  6,  and  6  ksi  \/In,  respectively. 

For  both  2024-T3  and  43140,  the  overload  interaction  zone  was  assumed 


to  be  equal  to  the  radius  of  the  plane  stress  plastic  zone, 


where 


(7) 


However,  a  and  hence  most  likely  depend  on  material  type  and  thickness 
Wei  et  (Ref  4)  presented  data  indicating  that  the  overload  interaction 

zone  in  Ti-6A1-4V  is  several  times  the  value  calculated  using  Eq  (7) . 
Accordingly,  the  predictions  for  Ti-6A1-4V  were  made  assuming  the  overload 
interaction  zone  to  be  approximately  four  times  the  plane  stress  plastic 
zone  radius  (  a  =  — ^  )  . 
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As  evident  from  Fig.  3,  the  Wheeler  model  as  modified  in  this  report 
can  predict  delay  cycles  following  a  single  overload  in  2024-T3,  4340, 
and  Ti-6AL-4V  within  essentially  a  factor  of  two  of  the  acttaal  experimental 
data.  For  all  cases  considered,  the  exponent  m  varied  from  1.07  to  4.34 
according  to  Eq  (6). 

The  model  at  present  appears  to  have  the  basic  elements  necessary 
for  a  generalized  crack  growth  analysis,  and  its  success  in  predicting  delay 
cycles  is  encouraging.  Further  refinements  are  being  made  toward  the  ultimate 
goal  of  predicting  crack  growth  under  a  general  load  spectrum. 
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TABLE  I 


Single  Overload  Delay  Data 


Material 

Reference 

R 

(ksijin) 

V 

(kslv'  in 

OL 

K 

max 

(Cycles) 

2024-T3 

3 

0.3 

6.S 

13.7 

2^01 

120,000 

9.9 

17.1 

1.73 

35,000 

20.4 

2.06 

87,500 

21.5 

2.17 

104,000 

12.2 

19.5 

1.60 

•  36,000 

23.7 

1.94 

107,000 

26.5 

2.17 

246,000 

14.4 

22.3 

1.55 

9,500 

30.6 

2.13 

398,800 

/  ■  ' 

32.4 

2.25 

278,500 

16.7 

27.8 

1.66 

21,000 

35.9 

2.15 

244,000 

5 

0.29 

12.6 

26.1 

2.07 

332,000 

21.0 

28.5 

1.36 

5,000  * 

36.0 

1.71 

52,200  * 

0.1 

8.2 

12.5 

1.52 

10,000 

14.9 

1.82 

32,000 

17.1 

2.09 

50,000  * 

*  Average 
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TABLE  I  (CONTINUED) 


Single  Overload  Delay  Data 


Material  Reference 

R 

K 

max _ 

(ksi  /in) 

OL 

(ksi,^) 

OL 

K 

max 

N  (Cycles) 

d 

2024-T3  5 

0 

13.7 

21.2 

1.55 

9,000 

28.5 

2.08 

72,000 

15.0 

30.0 

2.00 

98,000 

18.0 

27.0 

1.50 

13,000 

31.5 

1.75 

33,000 

36.0 

2.00 

113,000 

21.0 

31.5 

1.50 

14,000 

36.8 

1.75 

63,000 

42.0 

2.00 

258,000 

6 

0.05 

15.8 

23.2 

1.47 

5,000  * 

0.29 

21.0 

28.6 

1.36 

4,800  * 

0.67 

45.0 

52.4 

1.16 

5,400 

4340,  2 

0.1 

20.0 

40.0 

2.00 

35,000  * 

=120  ksi 

4340,  2 

0.1 

20.0 

40.0 

2.00 

4,000  * 

Oyg  =  220  ksi 

Ti-6A1-4V  4 

0 

14.0 

21.0 

1.50 

2,000 

28.0 

2.00 

7,500 

*  Average 
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TABLE  I  (C»NTINT]ED) 
Single  Overload  Delay  Data 


Material 

Reference 

R 

K 

max  _ 

(ksi 

^OL 

(ksi  \/ih) 

^OL 

Nj  (Cycles) 

^ax 

Ti-6A1-4V 

4 

0 

14.0 

31.3 

2.24 

15,000 

33.0 

2.36 

25,000 

33.6 

2.40 

30,500 

34.3 

2.45 

44,500 

7 

0.1 

18.5 

31.5 

1.70 

8,200 

37.0 

2.00 

11,200 

28.0 

42.0 

1.50 

2,400 

19.3 

38.6 

2.00 

10,100 

28.7 

43.1 

1.50 

2,500 
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TABLE  II 

CRACK  GROWTH  RATE  CONSTANTS 


da  ,  in/ eye  =  C  [AK,  ksi  '/in]^ 
dN 


Material 

^fer^nce 

R 

C 

n 

2024-T3 

3,5,6 

0.3 

1.48xl0"^ 

3.59 

2024-T3 

5,6 

0 

8.66xl0“^ 

2.68 

2024-T3 

6 

0.67 

-9 

3.21x10 

3.27 

4340,0  = 

ys  . 

120 

ksi 

2 

0.1 

-10 

1.38x10 

3.20 

4340. Oy3  = 

Ti-6A1-4V 

220 

ksi 

2 

0.1 

3.46x10“^° 

3.20 

4* 

0 

5.18x10”^ 

2.60 

Ti-6A1-4V 

7** 

0.1 

8.85x10“^ 

2.48 

*  AK  14  -  ksi  J  in 

**  AK  18  -  ksi  N/iiT 
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STRESS  INTENSITY  CRACK  LENGTH 


FIG.  1.  DELAY  CYCLES 


FIG.  2.  APPLIED  LOADING. 
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N^,  OBSERVED  DELAY 

CYCLES 

2024-T3,  Ref 

3 

□ 

2024-T3,  Ref 

5 

O 

2024-T3,  Ref 

6 

A 

4340,  Ref  2 

8 

Ti-6A1-4V,  Ref  4 

8 

Ti-6A1-4V,  Ref  7 

FIG.  3.  PREDICTED  DELAY 

CYCLES  VERSUS 

OBSERVED  DELAY  CYCLES 

